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The paper explores the interrelationship between the generalized Grad- 
Shafranov equation, or trans-field force balance equation, for steady MHD 
flows with an ignorable co-ordinate, and work by Imai on field-aligned MHD 
flows. The development of Imai, assumes at the outset that the fluid velocity 
V is parallel to the magnetic field B, and exploits an analogy with steady 
compressible irrotational flow in ordinary fluid dynamics. In Imai’s analysis the 
magnetic induction B is written in the form B = ob, where c oc 1/(M?, — 1), and 
M , is the appropriate Alfvén Mach number. Gauss’ law V.B = V. (ob) = 0 then 
plays à role analogous to the mass continuity equation in ordinary fluid 
dynamics, where o corresponds to the density of the pseudo-fluid. Imai’s 
analysis leads to a transonic equation for the field potential ¢ defined by 
b = V4. For a restricted class of flows the trans-field force balance equation 
formulation also leads to the transonic potential flow equation, but the 
assumption of an ignorable co-ordinate allows for the possibility of non-field- 
aligned flows with non-zero electric field potential ®,. The characteristics of the 
generalized Grad-Shafranov equation are related to the Mach cone and the 
group velocity surface for linear magnetosonie waves. The corresponding forms 
of the characteristics for the potential transonic flow equation in the (x, y) plane 
and in the (6,,6,) hodograph plane are discussed. Sample solutions of the 
potential transonic flow equation for radial, helical and spiral flows are obtained 
by means of the hodograph transformation, and are used to illustrate the 
differences between hyperbolic and elliptic flows. The potential transonic flow 
equation is obtained for the case of an ignorable co-ordinate z of a rectangular 
Cartesian co-ordinate system (x,y,z), and also for the case of flows with an 
ignorable co-ordinate $ of a spherical polar co-ordinate system (r,0,4). 
Astrophysical applications are briefly discussed. 


1. Introduction 


Models of steady MHD flows have been used extensively to model stellar 
winds and jets (see e.g. Heinemann & Olbert 1978; Tsinganos 1981; Sakurai 
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1985, 1990; Hu & Low 1989; Heyvaerts & Norman 1989; Camenzind 1986, 
1987), as well as models of bow shocks (see e.g. Spreiter, Summers & Alksne 
1966; Spreiter, Summers & Rizzi 1970; Spreiter & Rizzi 1974). Most models 
assume the existence of an ignorable co-ordinate, allowing the reduction of the 
MHD equations to the trans-field force balance equation (or generalized 
Grad-Shafranov equation), coupled with a sequence of integrals (e.g. the 
Bernoulli integral, the angular momentum integral, and the entropy integral 
for flows that have an ignorable azimuthal co-ordinate ¢ in spherical polar co- 
ordinates). Work by Imai (1960) and Seebass (1961) showed that, for the case 
of field-aligned flows, the MHD equations could be reduced to the potential 
transonic flow equation of compressible irrotational fluid dynamics by the 
introduction of the notion of an equivalent non-magnetized pseudo-fluid. The 
main purpose of the present paper is to explore the conditions under which the 
trans-field force balance equation for flows with an ignorable co-ordinate may 
be reduced to the potential transonic flow equation. The analysis leads to some 
interesting solutions of the MHD equations that are useful for delineating the 
differences between hyperbolic and elliptic flows, as well as the differences 
between sub-Alfvénic and super-Alfvénic flows. 

Section 2 introduces the steady MHD model employed. Section 3 considers 
the reduction of the MHD equations with an ignorable co-ordinate z of a 
rectangular Cartesian co-ordinate system to the trans-field force balance 
equation. A discussion of the characteristics of the trans-field equation reveals 
the conditions for the flow to be hyperbolic or elliptic. For hyperbolic flows the 
Mach angle between the characteristics and the projected flow direction on the 
(x,y) plane is determined, and the characteristics are related to the group 
velocity surface and Mach cone for linear magnetosonic waves. Section 4 shows 
the connection of the trans-field equation to the potential transonic flow 
equation for a class of flows for which the ‘source’ term in the trans-field force 
equation vanishes. The Mach number of the flow based on Imai’s pseudo-fluid 
analogy, M, is shown to coincide with the Mach number of the flow based on the 
analysis of the characteristics of the generalized Grad-Shafranov equation. The 
characteristics of the potential transonic flow equation in the (b,, bj) hodograph 
plane are discussed. For the case where gravity may be neglected the 
characteristics in the hodograph plane can be integrated explicitly, and used to 
determine the simple wave solutions. Section 5 uses the hodograph trans- 
formation to obtain sample solutions of the potential transonic flow equation 
for the cases of radial flow, helical flow, and spiral flows about the z axis. The 
analysis extends the work of Seebass (1961) to include the role of non-zero 
components of the velocity, V, and B,, not included in Seebass’ work. The 
representation of the flow and the characteristics in the hodograph plane are 
also discussed. In §6 the derivation of the trans-field equation for flows with an 
ignorable co-ordinate ¢ in spherical polar co-ordinates (r,6,¢) is carried out. 
The reduction of the trans-field equation formulation to the potential transonic 
flow equation is then developed for this geometry. Section 7 concludes with a 
summary and discussion. 
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2. The model 


The equations describing the steady inviscid flow of a high-electrical- 
conductivity MHD plasma consist of the dynamical equations for the gas, 
namely the mass continuity equation 


V.(pV) «0 (2.1) 
and the momentum equation 
pV.VV =IxB—pV®,—Vp, (2.2) 


supplemented by an appropriate equation of state relating the gas energy 
density to the entropy and density, and the Maxwell equations determining the 
inter-relationships between the electric field E, magnetic induction B and 
electric current density J. In the present analysis we assume that the first law 
of thermodynamics, or co-moving energy equation, is of the form 


V.V(p/p") = 0, (2.3) 


where y is the polytropic index of the gas. The choice of the polytropic index 
y in (2.3) may be used to simulate heat conduction if y is smaller than the 
canonical value of the specific heat ratio y,4. In (2.1)-(2.3) p, p and V denote 
the gas density, pressure and velocity respectively, and 4, denotes the 
gravitational potential. 

For a high-conductivity plasma E — — V x B, and Faraday's law in the MHD 
limit reduces to 


VxEzZ-Vx(VxB)-0. (2.4) 
The remaining Maxwell equations are Gauss' law 
V.B=0 (2.5) 
and Ampére's law 
J=VxB/n, (2.6) 


where 4 is the magnetic permeability. 
The Maxwell equations (2.4)-(2.6) yield Poynting’s theorem, or the 
electromagnetic energy equation in the form 


v5). - e (2.7) 


whereas the co-moving energy equation (2.3) and continuity equation (2.1) may 
be used to write the gas energy equation in the form 


v| vo J-v.v - 0. 2.8 
P G-p d a 
Taking the scalar product of the momentum equation (2.2) with V and adding 
the resultant equation to (2.7) and (2.8) yields the total energy equation in the 
fixed reference frame as 


v Lovan igs E - 0, (2.9) 


(y—-1)p 
where the terms from left to right represent the kinetic energy, enthalpy, 
gravitational and Poynting energy fluxes respectively. 


Exa) 
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3. Flows with an ignorable co-ordinate z 


We now briefly review the reduction of the above equations for the case 
where all physical quantities are independent of an ignorable co-ordinate z of a 
rectangular Cartesian co-ordinate system using an analysis similar to that of 
Tsinganos (1981). In §3.1 we provide a reduction of the steady MHD equations 
to the trans-field force balance equation or generalized Grad-Shafranov 
equation, followed by a discussion of the characteristics and Mach cone for 
hyperbolic flows. In §3.2 we discuss the connection between the group velocity 
surface for linear magnetosonic waves (see e.g. Whitham 1974), the charac- 
teristics, the Mach cone, and the differences between hyperbolic and elliptic 
flows. 

3.1. The trans-field equation 
Since z is an ignorable co-ordinate, the condition V.B = 0 implies that the 
magnetic induction B may be expressed in the form 
B= arem aser Bia) e, cB,-B,e, (3.1) 
Similarly, the mass continuity equation (2.1) implies that the fluid velocity V 
may be represented in the form 


.l(0p. Oy = 
= i (pe. ae) V(x, y)e, = V, + V,e,. (3.2) 
As in the development by Tsinganos (1981), the analysis leads to a 
generalization of the Grad-Shafranov equation for the magnetie potential 
A(x,y) that applies to steady flow situations. The magnetic field lines are 
described by the differential equations 


de dy _ d 
0A/0y  —ôðA/ðx B, 


z 


(3.3) 


These equations integrate to yield A(x, y) = const., for the projection of the field 
lines on the (x, y) plane. Similarly, the fluid streamlines have projections on the 
(x,y) plane described by the family of curves y(x, y) = const. 

Equation (2.4) implies that the electric field E may be written as the gradient 
of a potential (x, y) in the form 

E=—-VxB=-—V®,(z, y). (3.4) 

From (3.4), the z component of the electric field in the present model is zero. 
Using the representations (3.1) and (3.2) for V and B, one finds from the z 
component of (3.4) that E, = — (up) ! 0(w, A)/0(x, y) = 0, from which it follows 
that y = v(A). Also from (3.4), 


ôO, A) 
= 0, 3.5 
Q(x, y) em) 
which implies that ®, = 6,(4) is solely a function of A. Since y = (A), (3.1) 
and (3.2) imply that the components of V and B in the (x, y) plane are given by 


E.B = -V6,.B-— 


V 


p 


“A 
" rA )p,- “By (3.6) 
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where the subscript p denotes components perpendicular to the z axis, and 
a = y'(A)/u, The (x, y) plane Alfvén Mach number of the flow, M ,, is defined by 


2 
Mis. P BO um (3.7) 


(Vie Ba n 
where V ,, = B,/ (up)? is the (x,y) plane Alfvén velocity. 
The component of (3.4) parallel to VA yields 


V-5B. = (A). (3.8) 


The z component of the momentum equation (2.2) leads to the vanishing- 
Jacobian condition 0(B,—y’(A) V,,A)/0(z,y) = 0, which in turn implies the 
integral 
B 

yo——— 
* V'(A) 

where A(A) is an arbitrary function of A. 
Equations (3.8) and (3.9) can be solved for V, and B, to obtain 


— A(A), (3.9) 


De MAA- 9z) 
[LL A- OSA) 
Boma up lo (3.11) 


where M, is the Alfvén Mach number defined in (3.7). These equations clearly 
display a possible Alfvén singularity as M, > 1. Furthermore, if ®,(A) = 0 then 
V, = (a/p) B,, and the fluid velocity is then parallel to B, with V = (a/p) B. This 
latter case of field-aligned flow was considered by Imai (1960), whose analysis 
led to a potential transonic field equation similar to that occurring in steady 
irrotational fluid mechanics (see. e.g. Courant & Friedrichs 1948, chap. 4). 

Two further integrals can be obtained in a straightforward fashion. The first 
of these is the polytropic integral 


p = K(A)p", (3.12) 
which follows from integrating the co-moving energy equation (2.3), where 


K(A) is an arbitrary function of A. A further integral arises from the fixed-frame 
energy equation (2.9), namely 


©; (A) 
w+? eo, B, ET. = 
t ee T wa 
where £(A) is an arbitrary function of A. 
It is of interest to note that the energy integral (3.13) may also be cast into 
the form 


E(A), (3.13) 


E,(A) = E-M(A— 09) A] = MS (5-05) — e, (3.14) 
Equation (3.14) corresponds to the energy integral in a reference frame X’ 
moving with velocity V, = ok relative to the fixed frame. In the local frame 
X’ the flow is field-aligned and the electric field and Poynting flux vanish. 
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The component of the momentum equation (2.2) parallel to VA yields the 
generalized Grad-Shafranov equation or trans-field equation 


af ð (a.0A 6 {adA 1 /A @A 
5 lic ac) F5 a (at apt) 7 & iod 


Q = E'(4)- K(4) P — SA EY 
Y 


where 


LAN ATR AD) 


m UBI M- (3.16) 


Equation (3.15) may be further analysed by noting that p is an implicit 
function of A, VA, x and y, obtained by solving the energy integral (3.13) for p. 
The derivatives 0p/0x and 0p/0y on the left-hand side of (3.15) may be obtained 
by differentiating the energy integral (3.13). The net result of this procedure 
is a second-order quasi-linear nonlinear partial differential equation for A, of 
the form 


aA,,+bA,,+cA,,+d — 0, (3.17) 
where 
=A-V?,V2, b=—2V,, V,, Vi, 
a px’ p : pz py P I (3.18) 
c2 A-V3,V3, Am V-V, +a) + Vna] 


and a, — (yp/p and V, = B/(up): are the gas sound speed and Alfvén speed 
respectively. The explicit form of the coefficient d is unwieldy, and is given in 
the appendix. 

The characteristics of (3.17) are determined by solving the equation 


F = af? +b, é, +c% =0 (3.19) 
for £(z, y). Solving (3.19) for £,/, yields the equations 
Éy = Ss Es (3.20) 
where 
V, Voy V3 -- [A(V5 — A) 
—lmz "py "p p 
M AS (3.21) 


Thus the characteristics are real and (3.17) is hyperbolic provided that the 
discriminant 


ô = A(V5—A) = (Vi aj) (V5 — Vj) (V5 — Vi) (V5 — V2), (3.22) 


is positive. In (3.22) the speeds V, and V, correspond to the phase speed of fast 
and slow magnetosonic waves propagating along B,. The speeds V, and V, are 
obtained by solving the equation A = 0 for V,, namely 
Vi, = V5 +03 + [(V4 +03)? —4V4, a7]. (3.23) 
The speed V, in (3.22) is given by 
yV, = Vant = suh. (3.24) 
(Vata? (Ve+V5) 


For the case where B, = 0, V, = V,, and the speed V, then coincides with the 
group velocity of the slow magnetosonic mode propagating perpendicular to 
B (i.e. k.B = 0) at one of the slow-mode cusps. Since V, > V, > V,, 9 > 0 and 
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(3.17) is hyperbolic for V, in the ranges V, < V, < V, and V, > V;. The parameter 
ô <0, and (3.17) is elliptic for V,, in the ranges 0 < V, < V, and V, < Vp < V; 
The first-order partial differential equations (3.20) for £(x,y) have Cauchy 
characteristics 


iu; RR (3.25) 


The solutions of (3.25) for a given solution of (3.17) for A(x,y) define 
characteristic curves £(z, y) = const. and y(x, y) = const. in the (x, y) plane. 

To obtain further insight into the nature of the characteristics for the case of 
hyperbolic flows it is of interest to compute the Mach angle £ between the 
characteristic direction £ = V£/|V£| and the projected fluid velocity vector V, 
via the formula 

: IVE x V? 
sin? J = ————F.., 3.26 
vva pu 

To compute .Z, it is useful to note that the characteristic equation (3.19) may 
be cast in the form 


F = (A&—V$)|VE? + V3 VE x V? = 0. (3.27) 
From (3.26) and (3.27), we deduce that 


sin? of 2 (3.28) 
where we identify 
V: 
RA ie 
M? Wa (3.29) 


as the square of the effective Mach number of the flow. Equation (3.28) is 
analogous to the formula for the Mach cone angle s for supersonic flow past a 
stationary object emitting sound waves in ordinary gasdynamics (see e.g. 
Landau & Lifshitz 1987), in which the sound-wave disturbances are confined 
to the Mach cone. 

From (3.29), we find 


Vp- VP (Vp V5) 
(V4 aj) (V5 — Ve) 


For hyperbolic flow M? > 1, and the flow is ‘sonic’ (i.e. M = 1) when the flow 
speed V, matches either the fast or slow magnetosonic speeds V, or V,. 

An alternative form of the characteristic equation (3.19) may be obtained by 
converting to polar co-ordinates (r, ©), where (x, y) = r (cos ©, sin ©) gives the 
position vector in the (x, y) plane. In terms of r and O, (3.19) assumes the form 


V. V, A— VV; 
9E EFT POG = 0. (3.31) 


M:-1 (3.30) 


(A - V5 Vr) —2V 


This form of the characteristics is of interest in the consideration of the radial, 
helical and spiral flows considered in $5. 

In terms of the characteristic co-ordinates (£, y) obtained by solving (3.20) 
and (3.25), the generalized Grad-Shafranov equation (3.17) reduces to the 


canonical form 
PA, (E) A+) A, +d =0, (3.32) 
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where 
_ AE, E, A(A— V2) 
r= Ayy y. yi ; (3.33) 
ug) = a£ Obs Fee yy: (3.34) 


Equation (3.32) is the canonical form of the Grad-Shafranov equation in 
characteristic co-ordinates (é, 7) applicable for hyperbolic flows. 


3.2. The group velocity surface and Mach cone 


In this section our aim is to shed light on the nature of the characteristics (3.20) 
for hyperbolic flows, and the Mach cone described by (3.28) and (3.29), by a 
consideration of the group velocity surface for linear magnetosonic waves. We 
restrict the discussion to the case of field-aligned flows confined to the (x, y) 
plane. Related ideas on the characteristics for hyperbolic MHD flows are 
discussed by Sears (1960), Jeffrey & Taniuti (1964) and Polovin & Demutskii 
(1990) and references therein. We first discuss the group velocity surface for 
linear magnetosonie waves in a stationary medium using the approach of 
Whitham (1974) (see also Webb et al. 1993), and then go on to discuss the group 
velocity for magnetosonic waves in a steady flow and the Mach cone. 


3.2.1. The group velocity surface 
The propagation of linear magnetosonic waves in a stationary background 
medium is described by the wave equation 


o 2) 2* ys 20 oe 
s) =|F-(Wi+ay ov Valse =o. (3.35) 
where p for example represents the density perturbation (Lighthill 1960; 
Whitham 1974; McKenzie 1991). Here z, denotes distance along the uni- 
form background magnetic field direction. Plane waves with perturbation 
p cc exp [4(k. x — »t)] from (3.35) satisfy the magnetosonic dispersion equation 


o — (V4 t aj) wk? + V4 aj kz, k? = 0. (3.36) 
By using the identifications 


w=-S8S, k=VS, (3.37) 
where S is the wave phase, we obtain the wave eikonal equation 
St — (V4 +05) St (VS)? + V4 a2 82 (VS)? = 0. (3.38) 


This equation arises naturally in WKB expansions for MHD wave propagation, 
and plays a central role in defining the group velocity surface or eikonal surface 
for magnetosonic waves. 

Next we introduce cylindrical spatial co-ordinates about the magnetic field 
such that r = (x,,x,cos¢,2,sing), and co-ordinates k = k (cos ®, sin 9 cos b, 
sin sin) in k space (k = VS). The wave eikonal equation (3.38) has a 
complete integral in the form of a plane wave 


S = kx, cos +z, sin 9 — c(9) t] - S,, (3.39) 


where k, à and S, are regarded as arbitrary constants, and c(9) = w/k is the 
phase speed of the wave, satisfying the magnetosonic dispersion equation 


c! — (VÀ - a2) c? + V4 a2 cos? = 0. (3.40) 
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0? 30? 60° 90° 120° 150° 180° 


Ficure 1. Fast and slow magnetosonic wave speeds c, and c, versus wave normal angle à 
between k and B, for the case a = 1 and b — 0:9. 


Since w = c(#) k, the group velocity of the wave V, is given by 


V,- ce = (9) e, +e (8) ey, (3.41) 

where 
e, = cos 2 e, +sin 2 e,, (3.42) 
e; = —sin2e, t cos e,, (3.43) 


are unit vectors in the k and 2 directions respectively. 
The group velocity surface for the waves r — V,t can be constructed simply 
from (3.41)-(3.43) in the parametric form 


x, = Vat = [c(9) cos 9 —c'(9) sin 9] t| 
x = Viot = [c (9) cos Ò+ c(9) sin 9]t j 


by varying the wave normal angle 2. An alternative instructive way of 
constructing the group velocity surface described by (3.44) is to note that the 
wave eikonal equation (3.38) possesses a general integral as an envelope of plane 
waves for different # (Sneddon 1957). This general solution is obtained by 
determining 9 = P(x, z}, t) from the envelope equation S,;=0, and then 
substituting this value of 9 into the particular integral (3.39). Taking S, = 0, 
the simultaneous solution of the equations S = 0 and S, = 0 in the solution 
(3.39) yields the wave front (3.44) for a point source as an envelope of ‘plane’ 
waves. Thus the group velocity surface (3.44) may be thought of as an envelope 
of plane waves. 
Using the notation 
a=max{V,,a,}, bz min{V,,4a,}, (3.45) 


figure 1 shows the variation of the fast and slow magnetosonic wave speeds c; 
and c, obtained by solving (3.40) as 


e$, = Ha? +b? + [(a? +b°)? —40°b? cos? eP}, (3.46) 


(3.44) 


as functions of 9 for the case æ = 1 and b = 0.9. 
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FiaunE 2. Fast and slow magnetosonic wave eikonal or group velocity surfaces for a = 1, 
b — 09 and ( — 1. 


Figure 2 shows the group velocity surface described by (3.44) for fast and 
slow magnetosonic waves for a = 1 and b = 0:9. The group velocity is in the 
direction of the straight line vector OP = V, from the origin to a point P(2,, x2) 
on the wave eikonal surface. The wave vector k is normal to the surface. Taking 
S, = 0, the straight lines S = 0 described by (3.39) are tangent to the group 
velocity surface. 


3.9.2. The Mach cone 


To construct the Mach cone for supersonic flows, it suffices to determine the 
group velocity vector field i 
V,=V,+V (3.47) 
for magnetosonic waves at a fixed point in the flow, where V, denotes the group 
velocity of the waves in the fluid frame. In (3.47) V denotes the fluid velocity 
vector. In order to illustrate the general ideas, and to simplify the mathematics, 
we consider only the case of field-aligned MHD flow in the (x, y) plane. Hence 
the velocity field V = V, e, is taken parallel to the magnetic field B = B, in the 


(x,y) plane. 
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The Mach cone construction is illustrated in figure 3, which displays the 
group velocity surface f = V, in the fluid frame for the fast magnetosonic mode 
(we take t = 1 in the time-dependent formulae of §3.2.1). The group velocity of 
the waves in the fixed frame is obtained from the vector addition of v, and V 
via the Galilean transformation (3.47). Magnetosonic disturbances originating 
at the point Q in the flow are confined to the Mach cone P, QP,, where the 
limiting edge of the Mach cone P,Q is tangent to the group velocity surface. 

The first obvious geometrical aspect of the Mach cone from figure 3 is that the 
Mach cone angle £ and the wave normal angle 9 on the Mach cone are related 
b 

i +A =n, (3.48) 
a result that follows immediately from the fact that the wave vector k is 
perpendicular to V, on the Mach cone. 

From (3.39), the equation of the tangent P,Q to the Mach cone is described 
by the equation 

S = z, cos 9 - z, sin 9 —c(9) = 0, (3.49) 


obtained by setting S, — 0 and t= 1 in (3.39). Setting x, — 0 in (3.49) and 
solving the resulting equation for z,, we deduce that the vector QO representing 
the fluid velocity V, has magnitude 


c 
QOS) = i (3.50) 
Equations (3.48) and (3.50) now yield the basic formula 
1 
i -— 51 
sin S M (3.51) 
for the Mach cone angle «7, where 

M- E (3.52) 


is the Mach number of the flow based on the magnetosonic speed c(#) as 
measured on the Mach cone. The phase speed c(2) of the magnetosonie waves 
of interest on the Mach cone may be determined from the formula cos? = c/V, 
for the wave normal angle 2 from (3.50), and from the dispersion equation 
(3.40), to obtain 


2 
= Vy to} ane, (3.53) 
p 
Combining (3.52) and (3.53) yields the result 
yi 
M? =— 2 (3.54) 
VL—^ 


for the square of the effective Mach number of the flow, where A is defined in 
(3.18) with V, = V,,. Hence the analysis of the Mach cone based on the group 
velocity surface leads to the basic formulae (3.28) and (3.29) for the Mach cone 
angle % obtained previously from an analysis of the characteristics of the 
generalized Grad-Shafranov equation (3.17). 
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Figure 3. Construction of the Mach cone for the case of field-aligned flow in the (x, y) 
plane for hyperbolic flows with V, > V,. 


To determine the group velocity V, on the Mach cone, we note from (3.44) 
and (3.47) that 
V, = [ccos 9 —c'(9) sin 9 — V, Je, + [c (9) cos + csin 9]e,. (3.55) 
The unit vectors e, and e, are given by 
e= (A,, —A,,0)/B,, &= (4z, 4y, 0)/B,. (3.56) 
Noting that V, = (a/p)B, from (3.6), we obtain the formulae 
Vi, = —ALV,P— 18+ Vy], ] 


= 3.57 
Vi, = -AVM — ) x V] pen 
for the group velocity components on the Mach cone, where 
4 "m 3 1 
UA (3.58) 


Amp rcs i 
V (V5 +03) V5 —2V4 a2] 
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Ficure 4. Construction of the Mach cone for the case of sub-Alfvénic field-aligned flow in the 
(x,y) plane in which the fluid speed lies in the hyperbolic flow regime V, < V, < V,. 


Equations (3.57) may be used to show that the group velocity directions on the 
Mach cone coincide with the characteristic directions described by (3.20) and 
(3.21). A straightforward computation shows that 


ALT UE, (3.59) 


where £,/£, =¢, determines the characteristic directions for the generalized 
Grad-Shafranov equation (3.17). 

The Mach cone illustrated in figure 3 corresponds to a super-fast 
magnetosonic flow with V, > V, = max (V4,a,), in which the Mach cone occupies 
the region downstream of the source point Q of the waves. One can also 
construct a Mach cone for the slow magnetosonic flow regime V, < Vp < Vyp 
where V, = min (V,,a,), illustrated in figure 4. In this case the Mach cone P,QP, 
extends upstream of the wave source Q. It is not possible to construct a Mach 
cone in the elliptic flow regimes V, « V, « V, and 0 « V, « V, since the 
characteristics are then complex. 

A similar analysis of the Mach cone can clearly be carried out for flows not 
restricted to the (x,y) plane, but in this more general case the analysis is 
complicated by the fact that the waves of interest have their k vector restricted 
to the (x,y) plane, whereas B and V have components outside the (x, y) plane. 
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4. The potential transonic flow equation 


In this section we show the connection between Imai's analysis of field- 
aligned MHD flow and the generalized Grad-Shafranov equation (3.15). The 
potential transonic flow equation is derived in $4.1 for the case of flows with an 
ignorable co-ordinate z of a rectangular Cartesian co-ordinate system (z, y, 2). 
The characteristics of the potential transonic flow equation in both the (x, y) 
plane and the (b,, b,) hodograph plane are obtained in $4.2. 


4.1. The potential transonic flow equation 
Taking into account 


0 (x0A| | 0 («0A CA CA 
vxv,- [eas tay aay) lee Ber (Gat dale D 
the generalized Grad-Shafranov equation (3.15) may be written in the form 
pa B da 


V x [( M2) B].e, = npQ4- 7277. 


(4.2) 
An inspection of the form of Q in (3.16) reveals that Q = 0 for the special class 
of flows for which 

a —0, E'-AA'—-0, k’=0, A'—-05-0, (4.3) 
and (4.2) is then equivalent to 


Vxb=0, (4.4) 
where 
(1—43)B, 1B 1—M*, 
b = ——_* 2% =, g= ; 4.5 
(1—M 45) Boo T Bao 1—M?, VES) 


and M ,, and Bpo denote constant reference values of M , and B,. In terms of b, 
Gauss' law V.B — 0 assumes the form 


V. (cb) = 0. (4.6) 


Equations (4.4) and (4.6) are analogous to the equations of steady irrotational 
compressible fluid flow in which b plays the role of the fluid velocity and ø the 
role of the density. This analogy was first noted by Imai (1960), who considered 
the case of field-aligned flow. The present analysis indicates that this analogy 
also holds for cases including gravity (®, + 0) and for non-field-aligned flows, 
which occur in the presence of a non-zero electric field (57 + 0; see (3.10) and 
(3.11)). We note that the class of flows for which the conditions (4.3) hold 
includes the case where E, K, a, A and 6; are constant throughout the flow. 
Equations (3.10) and (3.11) also reveal that the choice A = ©; corresponds to 
a flow with B, —0, and V, — (4). This latter flow is free of the Alfvén 
singularity for V, and B,. More generally, (4.3) may be integrated to yield 


a=, E—1A-E,—-1A K=K,, A-05,-A,—95, (4.7) 
where quantities with subscript 0 on the right-hand sides of (4.7) are integration 


constants. If A, — 5, + 0 then V, and B, exhibit a singularity at the Alfvén 
point M, = 1. 
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For non-magnetized fluid flows, the energy conservation or Bernoulli integral 
(3.13) is of the form 


yrs f2+0,=8, (4.8) 


from which it follows that the adiabatic sound speed a, = (@p/ dp)? satisfies the 


equation 
2 nd 
a; = "(a ; (4.9) 


Hence for the pseudo-fluid described by (4.4) and (4.6) one may identify the 
square of the effective sound speed, a?, as 


a? = — ob (2) ] (4.10) 
Qc $; 
We note for later reference from (4.5) that 
dp p(M4—1) 
TE 4.11 
do oM} V 


To compute a’, we first note that the energy integral E from (3.13) may be 
written in the form 
_@ p wat(A-O2) Y y-1 / 19. 
B= Bet aga pt ty lake + Oot RATIOS. (4.12) 
In view of the conditions (4.3) imposed on E, a, A, K and 0; in the present 
analysis, differentiation of (4.12) yields 


m um = 1 ,M4-1, Vido db d 
o= (EAN) dd = Vi (re eJ ty (4.13) 


where M, = V,/a, is the sonic Mach number and a, = (yp/p)t the gas sound 
speed. The left-hand side of (4.13) is zero by the assumption E' — AA' = 0 in 
(4.3). From (4.10) and (4.13), we find 


Pa, nob rs. 44 
dogs- M'- yet (4.14) 


for the effective pseudo-fluid sound speed a and Mach number M, where A is 
defined in (3.18). From (4.14), the square of the pseudo-fluid Mach number 
M? = b*/a? is given by (3.29), which is just the square of the Mach number based 
on the characteristics (3.27) and Mach cone of the generalized Grad-Shafranov 
equation (3.17). Thus, from (3.30), the pseudo-fluid flow, like the Grad- 
Shafranov equation flow, is hyperbolic in the regions where M? > 1, namely 
V, > V, and V, < V, < V,, but is elliptic for V, in the ranges 0 < V, < V, and 
V, < V, < V;. The Mach number M may also be expressed in the form 


MM ERA 
TFAM E T T m 


dm v; 
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The continuity equation (4.6) for b can now be reduced to the transonic 
potential flow equation. First note that (4.6) may be written as 
«(s ôb ) b óc da 


zy p dr. 
F "a Y Oy 


x - 0. (4.16) 


However, from (4.13), 


ar [1i ea, nm 


Ox a\ * 9 yg V2 dx |’ 
20 = l Afe ien Aa) n9 on 
dy — MP dy  "Oy] Vz dy) 


Substitution of the results (4.17) into (4.16) results in the equation 


ab, ab db, Ob b? 
aop ee ei eot v4 22-5; 9. Ve, =0. (4.18) 


ax dy) Vi 
Since V x b = 0, b must also satisfy the equation 
Ob, ôb, 
LL. 4.1 
Ox dy aie) 


A natural way to incorporate the condition V x b — 0 is to write b as the 
gradient of a potential in the form b = V6, so that (4.18) reduces to 


2 2 
- $2) part- Gi) bry 725,9 6, 5, vé. ve 0. (420) 


Equation (4.20) is the usual transonic flow equation for non-magnetized 
compressible irrotational flow in the presence of gravity (see e.g. Courant & 
Friedrichs 1948, chap. 4; von Mises 1958, chaps 2 and 4). 

There is an extensive literature on the irrotational fluid equations (4.18) and 
(4.19) and the potential transonic flow equation (4.20) in the case where the 
gravitational force is negligible, since in this case (4.18) and (4.19) may be 
linearized by the hodograph transformation, by interchanging the role of the 
dependent and independent variables. Writing (b,, b,, 0) = (u, v, 0), one finds for 
$, = 0 that the hodograph transformation yields a linearization of the equations 
provided that the Jacobian of the transformation is non-zero: 


Q(u, v) 2 6(b,,, by) 
A(a,y) (a, y) 
For j — 0, one obtains the so-called simple-wave solutions. Some simple 
applications of (4.18) and (4.19) to ‘magnetic source’ and ‘magnetic vortex’ 
flow for field-aligned flow that uses the hodograph transformation and 


Chaplygin’s equation are given by Seebass (1961). 
Below, we discuss the characteristics of (4.18) and (4.19). 


j= (4.21) 


4.2. The characteristic equations 
The characteristics of (4.18) and (4.19) in the (x,y) plane in the hyperbolic 
regime consist of two families of curves C^ :a(x, y) = const., and C+: p(x, y) = 
const, satisfying the equation 


a? |dr? — |b x dr|? = 0, (4.22) 
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where dr = (dx, dy, 0) is the differential in the characteristic direction. Thus the 
curves C* and C satisfy the differential equations 


dy dx 
H^ DIM remm OA 
C+: 1 c. Ta’ (4.23) 
C: T Cs (4.24) 
where 
2. 42 
(, cito oer (4.25) 


2 542 
bi—a 


(Courant & Friedrichs 1948, chap. 2, $23, p. 47). From the characteristic 
equations in the form (4.22), the Mach angle o between the flow b and the 
characteristics is given by the equation 


Ibxdr?_ 1 


Wn 7 (4.26) 


sin? f = 


Thus, perhaps as might have been expected, the Mach angle for the pseudo-flow 
is just the same as for the Grad—Shafranov flow. The Mach angle  -> 0 and the 
characteristics are nearly aligned with the projected streamlines on the (x,y) 
plane as M — oo (i.e. for V, > V, and also as V, | V,). Similarly, as V, V, or as 
V,— V, one has M-i and ./-—90?, and the projected streamlines and 
characteristics are nearly orthogonal for weakly supersonic flows. 

The corresponding form of the characteristics '* and T in the (b,, b,) plane 
are 


db db b b.VO, dx 
Bo Im. f£ oILM g 
Ps da e da a*—bl V» da’ ise) 
-. db. b ^ b.V6,dx 
I: Bo tH ae Yid (4.28) 
An alternative form of the characteristic equations (4.27) and (4.28) is 
2 
.V® 
a? Kb — Ib. db — 720 qe qp 26; (4.29) 
p 


where db = (db,, db,, 0) and dr = (dx, dy, 0). 

One can regard the characteristic equations (4.23), (4.24) and (4.27), (4.28) as 
four partial differential equations for x, y, b, and b, as functions of a and f, 
provided that the transformation of co-ordinates from (x,y)- (a, f) is non- 
singular. The curves a(x, y) = const and f(x,y) = const in the (x,y) plane may 
then be identified with the characteristics C^ and C* respectively. 

In the absence of gravity (®, = 0), one can show that the characteristics C* 
and I", and the characteristics C^ and T* are mutually orthogonal. For ®, = 0, 
(4.18) and (4.19) are reducible in the sense that the characteristics I'* and T7, 
(4.27) and (4.28), then depend only on 6, and b,, and hence can be integrated. 
From the energy integral (4.12), b, a and M depend only on the density p, so that 
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b = b(p), a = a(p) and M = M(p). Inverting the relation b = b(p) then implies 
that a — a(b) and M = M(b) Using these facts and converting to polar co- 
ordinates in the (b,, b,) plane leads to the solution of the characteristics [+ and 
I" as 


d 1 db 
r+: 0—| WP- Ur — f — const, (4.30) 


z " 1 db 
D: 0+ (f*—1) 7, = a = const. (4.31) 


Equations (4.30) and (4.31) thus lead to the identification of the characteristic 
co-ordinates « and fl. I'* simple-wave solutions correspond to flows in which f 
is constant throughout the region of the (x, y) plane of interest. Similarly, in I7 
simple waves a is constant throughout the region of interest. From (4.29), the 
Mach angle £’ between the characteristics and b in the hodograph plane in the 
absence of gravity is given by 

_|b.db? 1 


cos? o^ = pede = Mi " (4.32) 
Hence £’ = 90+ .#° at points in the flow where the hodograph transformation 
is non-singular. 


5. Examples of magnetized flows with ignorable co-ordinate z 


In this section we provide sample solutions of the potential transonic 
equation (4.20) (with 9, — 0) by using the hodograph transformation as 
described by Courant & Friedrichs (1948), von Mises (1958) and Seebass (1961). 
The development generalizes the solutions obtained by Seebass (1961) to 
include the effects of non-zero components of the velocity and magnetic field, 
V, and B,, not included in Seebass’ work. For non-magnetized flows two 
fundamental solutions correspond to purely radial potential flow and circular 
flow. Other more complex flows are also possible, including the magnetized 
versions of spiral flows, Ringleb flows and simple waves. 

In $5.1 we briefly review some aspects of the hodograph transformation for 
(4.18) and (4.19) in the absence of gravity. In $5.2 we use these results to 
consider flows in which B, is purely radial in the (x, y) plane, of the form 


Lf 


B, = Bye, (5.1) 


p po 
where r denotes the radial distance from the z axis, and B,, and r, are constants. 
In $5.3 we consider the analogue of azimuthal circular flows and helical flows 


corresponding to the solution 
b- by es, (5.2) 
T 


for b. These two fundamental solutions are combined to yield examples of spiral 
flows in $5.4. 
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5.1. The hodograph transformation 


The essence of the hodograph transformation for the magnetized flow equations 
(4.18) and (4.19) (with ®, = 0) is that the equations may be linearized by 
interchanging the roles of the dependent and independent variables. Hence, in 
the present case, writing 


ehs v=b SUP (5.3) 


and assuming that the Jacobian j = O(u, v)/0(x, y) + 0, (4.19) and (4.18) may be 
rewritten as 
z,—y, = 0, (5.4) 


(a* —u*) y, (a? —*) x,  uv(z, y.) = 0, (5.5) 
where x = x(u, v) and y = y (u,v) are regarded as functions of u and v. The form 


of (5.4) suggests the introduction of a potential ®(u,v) such that 
x20, y-—6, (5.6) 


u? 


which ensures that (5.4) is automatically satisfied. Using (5.6) in (5.5) now leads 
to the linear equation 


(a®?—u?) ®,, + (à? — 7?) Duu -- 2uv6,, = 0, (5.7) 
vv 


for the potential ©. The corresponding equation for ¢ is the transonic potential 
equation (4.20) with ®, = 0, namely 


~ $3) Pza + (a? — 92) Pyy — 29. py Pay = 0. (5.8) 


Hence the hodograph transformation described by (5.3) and (5.6) leads to a 
linearization of the nonlinear transonic potential flow equation (5.8). Fur- 
thermore, the transformations (5.3) and (5.6) imply that ®(u, v) and ¢(z, y) are 


related by the Legendre transformation 
$ = ux+vy—¢. (5.9) 


In order to fully exploit the linearization provided by the hodograph 
transformation, it is useful to introduce the stream function y(x, y), where 


op 
ou , 0v ^ An (5.10) 
which automatically guarantees that V.B = 0. The stream function y(z, y) is 
clearly equivalent to the magnetic potential A(z, y) of (3.1). It also turns out to 
be useful to introduce the Legendre transform of y(x, y), namely W(w, v), where 


F(u, o) =r y $- = —ovr+ouy— Y, (5.11) 
and to develop corresponding equations for y and Y. 
To obtain solutions for purely radial flow and circulatory flow, it is useful to 
introduce polar co-ordinates (b, 0) in the (6,,6,) plane, where 


u =b, =bcosð, v=b, = bsinð. (5.12) 
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It then follows (see e.g. Courant & Friedrichs 1948) that ¢ and w are coupled 
by the equations 


ap _ bay 

06 o Ob’ (5.13) 
66 M*-10y 
b^ ob 6’ (5.14) 


where M = b/a is the pseudo-Mach number (4.14) for magnetized flows. The 
integrability constraint $,4— æ = 0 for (5.13) and (5.14) yields Chaplygin’s 
equation " dun i 

Oy 1-M'Oy 1-M Oy _ 

EE + bab + B ax 0. (5.15) 
The characteristics of Chaplygin’s equation (5.15) are given by the charac- 
teristics I', (4.30) and (4.31). It is also possible to obtain equations analogous 
to (5.13) and (5.14) linking ® and Y, namely 


or oo 
OV c(1—M9?)00 
7-5 28 (5.17) 
Elimination of ¥ between (5.16) and (5.17) yields the equation 
a? — b? 
aO + b? (5, + b®,) = 0, (5.18) 


for ®, which is the polar co-ordinate version of (5.7). Noting that ¢ = 6®,—9, 
(5.13), (5.14) and (5.18) yield the results 


by, = e(b0, —),, (5.19) 


= c (bo, + Dog), 5.20 
relating y and 6. Vo (b, + Pop) ( ) 
The Jacobian of the transformation from the (x, y) plane to the (u, v) plane is 


given by 


O(z, 1 ®,\" 
J= T 2 = | Poole t) (05) | (5.21) 
For the hodograph method to give a single-valued solution for w(x,y) and 
v(x, y), it is necessary to restrict the solutions of (5.8) to a region where J + 0. 
A line along which J = 0 is called a limiting line; and if a limiting line is crossed 
then the solution may become multi-valued. 


5.2. The radial B, solution 
The radial B, solution corresponds to the solution 
y n, (5.22) 


of Chaplygin's equation (5.15), where r, is a constant. Using the results (5.19) 
and (5.20) linking y and © leads to the corresponding solution 
> db 


$20(b)-n| Z (5.23) 
ob 
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for ®. Equations (5.6) then yield the equations 


= Tob, _ Toby 
ates, y= Oy (5.24) 


Upon noting B, = Bpo ab, (5.24) yield solutions for b and B, of the form 


x 


“Oe | (5.25) 
7 


which corresponds to the solution (5.1) for B,. 

The remaining aspects of the solution (5.1) for B, follow from the energy 
integral (3.14) (with ®, = 0), which may be written in the form of a potential- 
well equation as 


2Vp = E,—Utp), (5.26) 
where the potential U(p) is given by 
U(p) = KV,- epe ar (5.27) 


To obtain some idea of the physics implied by the potential-well equation 
(5.26), it is useful to introduce the dimensionless variables 
Ve, pel. quf. (5.28) 
Ojo Po r 
where a4, = (YPo / po. Py and r, denote the gas sound speed, density and radius 
at some canonical fixed point P, in the flow. The normalized form of the 
potential-well equation (5.26) is 


iV, = E,-U(p), (5.29) 
where 

Fu pr! V2 ,M3,(M3,— 1) 
U(p) = — + 2 ; 5.30 
a= (M3, — py n 
E, = Mio + + WV aco Mao, (5.31) 
My = », M= j^ > Vaso = Pen. (5.32) 

go Apo g0 


Taking into account the relationship (3.6) between V, and B,, and the solution 
(5.1) for B,, yields the expression 


V, = =e, (5.33) 


for V... Equations (5.29) and (5.33) implicitly determine the possible solutions 
for p — p(r), and the other physical quantities of interest then follow in a 
straightforward fashion from this relationship. The expression (5.31) for E, 
follows from the requirement that p = 1 at r — 1. 

In some of the examples we choose P, (i.e. r = 7,) to coincide with either the 
fast or slow magnetosonic point, so that the magnetosonic dispersion equation 


Voom (Vigo + ago) Vot Vipo aso =0 (5.34) 


6 PLA 52 
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Ficurs 5. The potential Ü(p) of (5.30) versus p for the parameters M ,, = 1:05, M,, = 2:5, 
Fizo = 069864 and y = 3. 


then requires that 
Lge = (Mo — 1) (Máy — 1) : (5.35) 
A0 
It is of course possible that the flow does not pass through one of the 
magnetosonic points, in which case (5.35) cannot be satisfied. In general, V, 
M ao and M,, are arbitrary parameters not satisfying (5.35). Hence, to obtain a 
solution of the normalized potential-well equation (5.27), it suffices to specify 
M30 Mso Vazo and y. The parameter 9; = 07/a,, needed to determine V,, (3.10), 
does not affect the potential-well equation (5.29), so that, without loss of 
physical insight, we set ©), = 0 in the calculations. 
An example of the potential well (5.29) is displayed in figure 5, where U(p) 
is plotted as a function of p for the case M, = 1:05, M,, = 2:5, V4.) = 069864, 
and y —$. From (5.35), this choice of parameters ensures that the point 
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4 5 6 


0 1 2 


3 
p 
Fiaurs 6. The square of the pseudo-Mach number M? versus p corresponding to the 
potential well of figure 5. For hyperbolic flows M? > 1, but for elliptic flows M? < 1. 


(p, T) 2 (1,1) corresponds to the fast magnetosonie point. For the above 
parameters, the energy constant E,, = 48941 is represented by the horizontal 
dashed line on the figure. The well naturally splits into a region of super- 
Alfvénic flow for which 0 < p < M*,, which is represented by the curve ABC in 
figure 5, and a region of sub-Alfvénic flow with p > M*,, which is represented 
by the curve DEFGH. The flow speed V, in the (z, y) plane is zero at the top of 
the well at the points C, D and H. From (5.26), (5.27) and (5.33), one can deduce 
that the points C, D and H correspond to V, = 0 and r = oo. The maximum flow 
speed V, in the super-Alfvénic region corresponds to the point A with p = 0 and 
r — oo. In the sub-Alfvénie region V, is à maximum at the bottom of the 
potential well at the point E. The corresponding plot of M? (where M is the 
effective Mach number) versus p is displayed in figure 6. Plots of p(7) and V, (7), 
for the super-Alfvénic flow region, obtained from (5.29), (5.30) and (5.33) are 


6-2 
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Ficurs 7. Plots of the density p and P, versus 7 for the two super-Alfvénic ‘radial’ flow 
solutions associated with the potential well in figure 5. There are two single-valued flows, 
namely the super-fast flow BA, and the near-Alfvénic flow branch BC, where the point B 
corresponds to the fast magnetosonic point with V, = V; 


displayed in figure 7. The corresponding profiles of p(7) and Vr) for the sub- 
Alfvénic flow region p > M%, are displayed in figure 8. 

To interpret figures 5-9, it is useful to compute the derivatives dV,,/dp and 
dV,/dr by differentiating the potential-well equations (5.29), (5.30) and (5.33). 
We obtain 


dV, _ (ag Va) (V5 Vp) 

LL NS 5.36 

dp PVV — Vian) Wen) 
AV, _ Vue Vaa) V3 - VIO T 
d r(V3— VUA VD | 


dr rV- VD VRV) 
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FiGURE 8. Plots of p and J, versus r for the sub-Alfvénic ‘radial’ flow solutions corresponding 
to the potential well in figure 5. There are two single-valued flows starting from the slow 
magnetosonic point F, where V, — V,, namely FED and FGH. 


where yog 
*2 — Ap%g 
P A al (5.39) 
In order to classify whether the flow is hyperbolic or elliptic, we recall from 
(3.30) that the Mach number M is given by 
(5 — V) Vp Vi) Vp ds 

S, ee 5.40 
QIEVDQL-VO' * Vita; MS 
For hyperbolic characteristics M? > 1, whereas for M? < 1 the flow is elliptic. 
The flow speeds V;, V,,, Vs, V5 and V, occurring in (5.36)-(5.40) are ordered 


according te V,» Vap > V$ > V, > Vy (5.41) 


The points on the potential well in figure 5 corresponding to the fast and slow 
magnetosonic speeds V, = V, and V, = V, are labelled B and F respectively; 


M=1+ 
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V, = V} at the point E, and V, = V, at the point G. From figure 6, the flow is 
hyperbolic in the super-fast magnetosonic region V, V, V,,, which 
corresponds to the curve segment AB in figure 5. For V,, < V, < V, (curve 
segment BC) the flow is elliptic. In the sub-Alfvénic flow regime the flow is 
hyperbolic for V, « V, « V,, corresponding to the segment FG in figure 5, and 
is elliptic on the segments DEF (V4, < Vp < V?) and GH (0 < V, < V,). On the 
plot of M? versus p in figure 6, M? > oo as p— p, = 42333, where V, = V,. Note 
from (5.40) that M? < 0 in the region p > p,, and that M? = 1 at the fast and 
slow magnetosonic points B and F in figure 6. 

The plots of p and V, versus 7 for the super-Alfvénic flows in figure 7 are 
double-valued, and the derivatives dV,/dr and dp/dr diverge at the fast 
magnetosonie point B, where V, = V, (see also (5.37) and (5.38)). Hence, to 
obtain a single-valued solution, it is necessary to restrict the flow to be either 
supersonic in the fast magnetosonic sense with V, > V, (segment AB), or 
subsonic with V,, < Vp < V, (segment BC). One can also see the possibility of 
multi-valued solutions from the Jacobian (5.21) of the hodograph trans- 
formation. Using the solution (5.23) for ® in (5.21) yields the expression 


Quy) uf Y OUS DP - 1) 
enG) TEZI” (5:32 


for the Jacobian of the hodògraph transformation. Hence J = 0, and limiting 
curves occur in the (x, y) plane at the radii r = 7, and r = r, at the fast and slow 
magnetosonic points, where V, = V; or V, = V,, and M? = 1. The Jacobian (5.42) 
diverges when V, = V, or when V, = V,,. From (5.33), the constant radial mass 
flux for the flows in figures 7 and 8 is given by 


m = 2nrp Vp = 2mr, Po Vs. (5.43) 


Thus the supersonic (V, > V;) and subsonic (V, < V;) flows in figure 7 have the 
same mass flux and energy constant Em. The supersonic branch consists of a 
high-speed low-density flow, whereas the subsonic branch consists of a low- 
speed high-density flow. 

Figure 8 illustrates the radial profiles of p(r) and V,(r) for the sub-Alfvénic 
flows FED and FGH. There are two possible single-valued solutions starting 
from the slow magnetosonic point F, namely the low-density branch FED with 
Vi > V, > V, and the high-density branch FGH with 0 < V, < V,. The branch 
FED is supersonic in the slow magneto-acoustic sense, but is an elliptic flow 
with complex characteristics. Starting from the slow magnetosonic point F, V, 
first increases until it achieves its maximum at the point E, where V, = Vj, and 
then decelerates with increasing r. The main source of energy driving the flow 
on the segment FE is the gas enthalpy, with the density changing from 
p = p, = 3:828 at F to p = p, = 1:2857 at E. At larger radii r > rą, beyond E, the 
radial flow decelerates and V, increases at the expense of V,, (see figure 9), and 
the enthalpy remains roughly constant. The subsonic, slow magneto-acoustic 
branch FGH with V, < V, consists of a hyperbolic section FG with V, < Vp < V, 
followed by an elliptic flow regime GH at larger radii with 0 « V, « V,. The 
kinetic flow energy 3V;, associated with the radial flow is converted to enthalpy 
yp/[(y—1)p] as the flow moves outward, with a relatively minor role being 
played by the z component of the kinetic energy 3/2. 

Further features of these solutions are presented in figure 9, which displays 
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Ficure 9. Radial profiles of V,, B, p, V, B,, M, and M, for ‘radial’ flow solutions 
corresponding to the potential well of figure 5. 


the radial variation of V,, V,, B, B, and p (note that B = B/B,, and 
V = V/a,,). Also shown are the radial variation of the fast magnetosonic Mach 
number M, = V,/V; for the super-Alfvénic flow regime and the variation of the 
slow-mode Mach number M, = V,/V, for the sub-Alfvénic flow regime. The 
normalized forms of V, and B, from (3.10) and (3.11) are 


coc om ca Malis) 

V, = Drt Vu, 5.44) 
E Azo Mh.—p ( 
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FicurE 10. The potential well Ü(p) and the square of the pseudo-Mach number M? versus p 
for flows with V, = B, = 0, in which M, = 1, M, = 25, Vizo = 0 and y =$. 


At sufficiently large radii r > r; (r = rj at B where M, = 1) the radial flow speed 
profile V,(r) for the supersonic slow mode branch FED is almost identical with 
the subsonic fast-mode branch BC. Both the flows BC and ED are almost 
Alfvénic in the sense that M, %1 for these flows. For these solutions V, 
decreases and |V,| increases with increasing r. However, V, > 0 and B, > 0 for 
the super-Alfvénie solution BC, whereas V, « 0 and B, « 0 for the sub-Alfvénic 
solution FED. The density profiles p(r) show that the enthalpy and density 
variations play a more dominant role in the sluggish sub-Alfvénic flows than in 
the lower-density super-Alfvénic flows. 

A further example of the potential-well solution (5.29) with $5 = V, = 0, 
M,,=1, M,,—25 and y —$ is displayed in figures 10 and 11. Since 
r = Vi.) = 0, V, = B, = 0. This example corresponds to the type of solution 
considered by Seebass (1961) in which the fluid velocity is restricted to the (x, y) 
plane. Figure 10 (a) shows the potential-well plot of U(p) versus p. Unlike the 
potential well in figure 5, there is no singularity at the Alfvén point p = M*, = 1. 
Since V,, = 0, the fast magnetosonic point occurs when V, = V,, = V, at B. The 
slow magnetosonic point occurs at the point C, where V, = V, = a,. At the point 
D the flow speed V, matches the slow magnetosonic cusp speed (i.e. V, = V,). 
The solution is degenerate in the sense that two of the characteristic MHD 
speeds, namely the Alfvén and fast magnetosonic speeds, coincide in this 
example. From the plot of M? versus p in figure 10 (b), the solution is hyperbolic 
in the super-Alfvénic flow regime p < M*,. In the sub-Alfvénic flow regime the 
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FrcunE 11. Radial profiles of p and V, for the ‘radial’ flow solutions corresponding to the 
potential well in figure 10. There are two single-valued solutions starting from the gas sonic 
point C, where V, = a,, namely CBA and CDE. 


solution is elliptic on the segment BC, where a, < V, < Vap, hyperbolic on CD, 
where V, < V, < a,, and elliptic on DE, where 0 < V, < V,- The plots of p(7) and 
V,(7) in figure 11 show that the solutions join in a continuous fashion at the gas 
sonic point C. Formally, as V4,, ^ 0, one has V5 — V4,, and (5.37) and (5.38) 
reduce to 


ONG oc Pig. (5.46) 


A eec P pr (5.47) 


Hence |dV,/dr| oo and |dp/dr|- oo as V,—a,. Thus there are two possible 
single-valued solutions in figure 11, namely the supersonic branch CBA with 
V, > a, and the subsonic branch CDE with V, < a,. The solution CBA consists 
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of an elliptic section CB followed by a hyperbolic flow regime BA. Similarly, the 
subsonic solution CDE consists of a hyperbolic flow region CD followed by an 
elliptic flow region DE at larger r. The near-Alfvénic flow solutions ED and BC 
in figure 9 do not exist for V,,) = 0. 


5.3. Helical flows 


In this section we consider the character of helical and circular flows associated 
with the solution (5.2) for b. This corresponds to the solution 


D =b 7,9 (5.48) 


of the linear equation (5.18) for ®. Equations (5.6) then yield the solution (5.2) 
for b in the form 
b, = boty, by m — bre, (5.49) 
r r 
which corresponds to clockwise rotation in the (x, y) plane of the vector field b 


if bare > 0. The corresponding solution for the stream function y may be 
determined from (5.19) and (5.20) as 


y= -rr fo, (5.50) 


where c is the pseudo-fluid density (4.5). 
Taking b, = — 1, (5.2) and (4.5) yield the solutions 


Mairs 
Po yir ® 


for b and B,, where e, is the azimuthal unit vector in the (x,y) plane. The 
remaining aspects of the solution then follow from the potential-well equations 
(5.26)-(5.32) used in the discussions of the purely radial flow solution in $5.2. 
Since the normalized magnetic field component B, = B, in the (x, y) plane is 
related to the corresponding velocity component f = = V, by the equation 


b- Teg, B,=B (5.51) 


B,=B, =h (5.52) 


(see (3.6)), we may identify the behaviour of the azimuthal field with that of the 
azimuthal mass flux. Equations (5.51) and (5.52) and the potential-well 
equation (5.29) implicitly yield the dependence of the density p on the 
cylindrical radius r via the equation 

a2 MO — 1)? 

ga, - py En- OP id 
The components of V and B along the z axis follow from (5.44) and (5.45). 

Figure 12 illustrates an example of a helical-flow solution for the case 

M= 1:1, My 220, y 2$ and V,,, = 072157, where the value of V,,, is 
determined by the normalization equation (5.35). The figure shows plots of the 
velocity field components V, and V,, the magnetic field components B, and B,, 
the density p, and the fast and slow magnetosonie Mach numbers M, and M, as 
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Ficure 12. Radial profiles of P, Bs, p, V,, B., M, and M, for helical flow solutions with 


2 


M, = 1-1, Mo = 20, Vi; = 072157 and y = $. 


functions of 7. The points A, B and C on the super-Alfvénic flow solution 
(V, > V4,) and the points D, E, F, G and H on the sub-Alfvénic flow branch 
have the same physical meaning as in the potential-well diagram illustrated in 
figure 5. Hence C, D and H signify the points at r = oo, where V, — 0, and 
correspond to the top of the potential well. The points B and F correspond to 
points in the flow where the azimuthal flow velocity V; matches the fast and 
slow magnetosonic speeds V, and V, respectively. On the sub-Alfvénic branch 
V, = V} at E and V, = V, at G. As in the radial flow solution examples, the flow 
is hyperbolic in the super-Alfvénic flow regime AB, where V, > V;, and for the 
slow-mode regime FG, where V, < Vp < V, but is otherwise elliptic. 
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To obtain some physical insight into the character of the solutions, it is useful 
to note that the momentum equation (2.2) for this class of flows reduces to the 
radial momentum equation 


B- -2-$0525)-5 


= 54 
n Je (5.54) 


P, dr dr 


Hence the centripetal force on the fluid, —pV2/r, is balanced by the gas and 
magnetic pressure gradients and the magnetic tension force, — B2/ur. The 
balance between the different force components in (5.54) is significantly 
different in the different flow regimes. Thus, for example, in the elliptic flow 
regime 0 « V, « V, (GH) the gas pressure gradient, —dp/dr, and the magnetic 
tension forces, — B2/ur, are inwards, but the magnetic pressure gradient force, 
— d[(B2 - B2)/2,]/dr, is outwards. On the other hand, for the supersonic fast- 
mode flow regime V, > V, (AB) all the forces on the fluid element are inwards. 

Further aspects of the solution may be assessed from the radial derivatives 


dB, B,(V5—Vp(V5— V) 


dr o7 A Ea (YS um 


dp p VyV?— Vig) _ 
dr rV +) V-V (9:86) 


2 —yx*2 
dV, _ _ V,(ag Viz) (V5 — V5") (5.57) 


ERU ET TN 


Hence B, = B, and the azimuthal mass flux pV, (see (5.52)) have maxima when 
the azimuthal flow speed V, matches either the fast or slow magnetosonic 
speeds V, or V, at the points B and F; V, is maximum in the slow-mode regime 
at the point E, where V, = V?. Each of the derivatives in (5.55)-(5.57) diverges 
at the point G, where V, = V, 

The maximum azimuthal speed for these flows occurs on the fast-mode 
branch at the point A, where cavitation occurs (i.e. p— 0), at a radius r = fmin» 
which is the value of r at p — 0 given by (5.53). The flow is restricted to the 
region r> Tmin: The V, velocity component at r= fmin is non-zero, and is 
determined by setting p = 0 in (5.44) (note that ©; = 0 is assumed in figure 
12). The azimuthal flow speed V, decreases monotonically and V, increases with 
increasing radius on the fast-mode branch ABC. The velocity V > V, e, as r> oo, 
so that the flow is essentially non-rotating at large radii. 


5.4. Spiral flows 


Solutions for spiral flows can be obtained by linear superposition of the 
potentials ® corresponding to radial, (5.23), and circular, (5.48), flows. The 
resulting potential ® is of the form 


o6-—n (cosa) [in 9|. (5.58) 


where the first term on the right-hand side represents purely radial flow, and the 
0-dependent term represents circular flow. Since (5.18) for ® is linear, a linear 
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superposition of independent solutions is again a solution. The corresponding 
stream function y is readily determined from (5.19) and (5.20) as 


v=, [csin Qo) IE: + (eos a) J : (5.59) 
From (5.6) and (5.58), we find 
r= 7a sin a) by + (cos Zt 
— To ; by 
7g — (sin a) b, + (cos Oy) s (5.60) 
Equations (5.60) are readily inverted to yield 
T X 5. 
b, = “(cos %) —— (sin Xo) j| 


b, = “alsin a) + (cosa), (5.61) 


where 7? = z?-Fy?. Alternatively, using polar co-ordinates (x,y) = r(cos®, 
sin ©), (5.61) give the solution 


sin a 
b= e, eo (5.62) 


for b, where F = r/r, is the normalized form of r. 
Having established the solution (5.62) for b, (4.5) and (5.45) give the solution 
for B = B/B,, in the form 


=- B 5; Q-1p. Vaio MaMo 1) 
B= — ro +B A0 + Az0^^ A0 AO : 5.63 
T a "e (M5.,—p)r is M5, — p) 2 ) 
where EN 
(Bro Boo) = (cos o, sin œo) (5.64) 


are the normalized components of B, at the canonical radius r= ry. The 
corresponding normalized velocity field V = V/a,, for the case ®p =0 is 
parallel to B, and is given by 


V =M. (5.65) 


As in §5.2, we choose the radius r = r, to coincide with the fast magneto- 
sonic point, so that V, is determined from (5.35). The polar co-ordinate form 
of the transformations (5.61), with b = b(cosO,sin0) yields the equivalent 


transformations 
© = 0—tan ! (c tana), (5.66) 


cos? a, 
o’. 


P 


7b? = sin? a + (5.67) 


Since ø = 1 at r = r,, (5.66) yields 
Ay = 0,— 9,, (5.68) 
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Ficure 13. Radial profiles of V,, B, and p for spiral flow solutions with M ,, = 1:1, M, = 2:0, 
Vizo = 012157 and y = $. The profiles have vertical tangents at the limit lines at the points 
B and H, where the radial flow speed V, — c, with c being the appropriate magneto-acoustic 
phase speed on the Mach cone. 


so that a, may be identified as the angle between B, and the radial direction at 
Foz dy. 

The remaining equation needed to determine the solution is the energy 
equation (5.29). From (5.29) and (5.65), the dependence of the density p on 
radius 7 is determined implicitly from the potential-well equation in the form 


p Mi [B5 + Be (Mao— 1)? £/00.,— py] 


d 2 5.69 
op En - O(p)] (5:69) 
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The resulting solutions for p, V and B depend only on radial distance 7 from the 
z axis. 
The streamlines of the flow, 


dr rdQ dz 
using the solution (5.65) for V may be integrated to yield 
" dr 
9— 6, = (tan a) ges (5.71) 
To 
P M 
ZZ, = VazoM ao oF dF. (5.72) 


M, COS a Jy. 


The projection of the streamlines (5.71) on the (x,y) plane is of course 
equivalent to the equation w = const, where y is the stream function (5.59) — a 
result that may be verified from the transformations (5.66) and (5.67). 
Examples of the spiral flow solution (5.62)-(5.69) for the parameter values 
Ma = 1:1, Myy 2, Vazo = 072157, and a, = 45? are given in figures 13 and 14. 
Figure 13 shows the dependence of V,, B, and p on radial distance 7 from the 
z axis, whereas figure 14 shows the variations of the magnetic and velocity 
vector field components (B,, Bo, Bz) and (V,, Vo, V,). To interpret figures 13 and 
14, it is useful to determine the derivatives dp/dr, dV,/dr and dB,/dr from 


(5.63)-(5.65) as 
dp oV- VV: 


CL p o app 
T - 1 (5.73) 
d (a2 A )( 1) 
p— png a vp p-. 
u HI : l i ; (5.74) 
dB B,(V2,— V3 (V2 — V2) 
p. —»'np 'f/NVp "si 
T : ; (5.75) 
where D -A-Và ys, Ah (V2, — V3) (Và — V3). (5.76) 


In (5.73)-(5.76), V, and V, denote the fast and slow magnetosonic speeds (3.23), 
and V7 is defined in (5.39). 

To obtain some physical insight into the singularity D = 0, in (5.73)-(5.75) it 
is useful to note from the Mach cone analysis of $3, that the Mach number M 
of the flow from (3.29) and (3.52) is given by the formulae 


y? V3 
M=2=—_ 2, (5.77) 
è y$-a4 


where c denotes the appropriate magnetosonic phase speed measured on the 
Mach cone. From (5.76) and (5.77), we find 

A= P-e), D = VV e). (5.78) 
Hence the radial derivative |dp/dr|-> oo and D0 when the radial flow speed 
V, matches the magnetosonic phase speed c on the Mach cone. At this point in 
the flow the characteristic equations (3.31) reduce to 


ter | =0. (5.79) 
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Fieure 14. Radial profiles of the velocity field components (V,, Vo, V.) and magnetic field 
components (B,, Bo, B,) for spiral flow solutions with the same parameters as figure 13. 
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One of the solutions of (5.79), namely £, = 0, shows that locally at the radial 
sonic point, where D = 0 and V? = c?, this characteristic is purely azimuthal (i.e. 
E = O at this point) and the k vector of the magnetosonic wave of interest is 
purely radial. In terms of the pseudo-fluid velocity b and sound speed a, one ean 
obtain equivalent formulae for 4 and D, namely 


_ Vpb) p Vba) 


b? b? 


— 


A D (5.80) 


~ 


Hence at the radial sonic point D = 0 and b? = a’. 
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Ficure 15. Examples of streamline projections on (a) the (x, y) plane and (b) the (z, r) plane 
corresponding to the super-Alfvénic spiral flow solution BA of figure 13. (c) shows the 
characteristics I and T™ in the (b,, b) hodograph plane, and the streamline S that meets the 
characteristics on the limit circle at b — b,. (d) shows the corresponding characteristics C* and 
C- and streamline S in the (x,y) plane. 


It is evident from figure 13 that at the points B and H, where V? — c? and 
D = 06, the fast and slow magnetosonic solutions become double-valued. This 
fact is related to the vanishing of the Jacobian (5.21): 

O(x,y) | r*(bi—a?) r(A — V? V2) 
Se 5.81 
00,6) a PTAA V- VA um 
at the radial sonic point, corresponding to a limiting line in the hodograph 
transformation. The Jacobian (5.81) diverges as V,>V,, where c>0 and 
WP ^ oo. 


5.4.1. Streamlines and characteristics for spiral flows 


Although the results of figure 14 contain the basic information on the 
variation of V and B for the spiral flow, a more vivid characterization of the 
flow is obtained from the flow streamlines (5.71) and (5.72) and from the 
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b, x 
(c) (d) 
Figure 16. Same as figure 15, except for the near-Alfvénic spiral flow branch BCD of figure 


13. The circles b = b, and b = b, correspond to the limit circle and fast magnetosonic point 
respectively. 


hodograph plane characteristics [+ and I^ described by (4.30) and (4.31). The 
streamlines (5.71) and (5.72) and the characteristics (4.30) and (4.31) are most 
simply evaluated by converting the integrals involved to integrals over the 
density p. 

The streamlines and characteristics for the fast magnetosonic spiral flow BA 
of figure 13 (obtained by numerical integration) are displayed in figure 15. 
Figures 15(a) and (b) show the streamline projections on the (x,y) and (z,7) 
planes. Figure 15(a) shows that the spiral flow branch BA is almost a radial 
flow, with the circle at radius r = r, representing the limit line (i.e. the point B 
in figure 13). 

Examples of the characteristics I'* and T in the (b,, 6,) hodograph plane are 
shown in figure 15(c). The circles at b = b, and 6 = b, correspond to the limit 
line r = r, and r = oo respectively. The dashed curve S represents a fluid stream- 
line in the hodograph plane. The constants in the integrals (5.71), (5.72) and 
(4.30), (4.31) have been chosen so that the streamline S and the characteristics 
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Figure 17. A unified view of the characteristics and streamlines S of figures 15 and 16 in the 
hodograph plane. The flow BA lies outside the limit circle b = b,, whereas the flow BCD lies 
inside the limit circle. 


meet at the limit circle b = b,. By varying the integration constants, the hyper- 
bolic flow region b, < b < b,, can be covered by a net of characteristics I'* and 
I^. The behaviour of the flow streamlines and characteristics is complicated in 
the immediate vicinity of the limit line, where the hodograph transformation is 
singular. A discussion of the character of the hodograph map and characteristics 
in the vicinity of the limit line for the case of pure gas dynamics is given by 
Courant & Friedrichs (1948, chap. 2, $30). The characteristics [ and T" in 
figure 15(c) are visually similar to the epicycloid characteristics for gasdynamics 
discussed by Courant & Friedrichs (1948, chap. 4). The corresponding 
characteristics C* and C^ and streamline S in the (x,y) plane are displayed in 
figure 15(d). Note that the characteristic C* is tangent to the limit circle at 
r=). 

The second possible solution for super-Alfvénic flow in figure 13 is the 
solution branch BCD. The point B corresponds to the limit line where V, = c, 
and the point C represents the fast magnetosonic point where V, = V;. The 
projections of the streamlines on the (x, y) and (z,r) planes for this spiral flow 
are shown in figures 16(a) and (b). Unlike the flow in figure 15, this near- 
Alfvénic flow has a pronounced spiral structure. The flow is hyperbolic between 
the two circles r = r, (the limit line) and r = r; (the fast magnetosonic point). 
Figure 16(c) shows the characteristics F* and P^ and the streamline S, which 
meet at the limit circle b= b, The inner circle b =b, represents the fast 
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FIGURE 18. Streamlines and characteristics for the sub-Alfvénic spiral flow HGFE of 
figure 13. The flow is hyperbolic for b, < b < b, in (c). 


magnetosonic point. The characteristics C^ and C^ in the (x, y) plane are shown 
in figure 16 (d), in which the characteristic C* is tangent to the limit circle r = r,. 
A unified view of the two super-Alfvénic spiral flow solutions is presented in 
figure 17, which shows the characteristics [* and I" and the spiral flow 
streamlines S in the hodograph plane for the flows BA and BCD depicted in 
figures 15 and 16. The points A and D correspond to r = oo, and the two single- 
valued flows extend from the limit line b = b,. 

Figures 18 and 19 show respectively the sub-Alfvénic solutions HGFE and 
HIJ in figure 13. The main point to note is that the spirals rotate in the opposite 
sense to the super-Alfvénic spirals in figures 15 and 16, and V, and B, are along 
the negative z axis. In order to resolve the structure of the streamlines and 
characteristics, figure 20 shows an expanded view of the hyperbolic flow region 
in the hodograph plane. The curves labelled b = b,, b = b, and b = b, represent 
the slow magnetosonic circle (V, = V,), the limit circle (V, = c) and the slow 
magnetosonie cusp circle (V, = V,) respectively. The streamlines 8, and 5$, 
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FIGURE 19. Same as figure 18, except for the sub-Alfvénic spiral flow solution HIJ of 
figure 13. The flow is hyperbolic for b in the range b, « b « b, in figure (c). 


associated with the flow in figure 18 are tangent to the characteristic I'* at the 
limit circle. 


6. Flows with an ignorable azimuthal co-ordinate ¢ 


In this section we indicate the extension of the results of §§3 and 4 to the case 
of steady MHD flows with an ignorable co-ordinate ¢ of a spherical polar co- 
ordinate system (r,0,9). The corresponding version of the generalized 
Grad-Shafranov equation, or trans-field equation, analogous to (3.15) has been 
used extensively in modelling magnetized wind flows and jets (see e.g. 
Heinemann & Olbert 1978; Tsinganos 1981; Sakurai 1985, 1990; Camenzind 
1986, 1987; Hu & Low 1989; Heyvaerts & Norman 1989). We first outline the 
derivation of the trans-field equation for the above geometry, and then 
determine under what conditions the trans-field equation coupled with Gauss’ 
law V.B = 0, leads to the potential transonic flow equation. We remark that the 
conditions for the potential transonic flow to apply are sufficiently stringent 


182 G. M. Webb, M. Brio and G. P. Zank 


4.5 4-6 4-7 4.8 

b, 
Ficure 20. Expanded view of the characteristics I'* and I for the spiral flows of figures 18(c) 
and 19 (c). The curves b = b, b = b, and b = b, are segments of circles in the hodograph plane 
corresponding to the slow magnetosonic point, the slow-mode limit circle and the 
magnetosonic cusp point where V, = V.. S, and S, are the streamlines of figures 18(c) and 
19(c), which both issue forth from the limit line. - 


that it is not possible to satisfy the usual regularity conditions applied in wind 
theory at the Alfvén critical surface. Hence this approach will not be of any use 
for models of winds near the Alfvén surface, but could be of use for sub-Alfvénic 
siphon flows in the solar atmosphere, or for subsonic and supersonic flows past 
a stationary object (e.g. the magnetosphere of a planet). Solutions of the 
potential transonic equation should also be of use in describing the collimation 
of radio jets by the Lorentz force. 


6.1. The trans-field equation 


Since ¢ is an ignorable co-ordinate, Gauss' law V.B — 0 and the continuity 
equation (2.1) imply that B and V may be expressed in the form 


A 
B=Vx (ae t B,(r,0)e, = B, B,e,, (6.1) 
vedo v es E (6.2) 
Tao’ Araing 4] * 659 7 Vat "aee 


The projection of the magnetic field on the (r, 0) plane is given by the contours 
of A(r,@) = const, whereas y(r,0) = const gives the projection of the fluid 
streamlines. It is of interest to note that B, = V x (AV¢) = VA x Vó. 
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The ¢ component of Faraday’s law in the form E = —V x B = —V®,(r, 6) 
gives 

= 1 Oy, A) _ 
$ ppr'sin*Ó Q(r,O) — 0, (6.2) 
and hence y = (A). Similarly, since E.B = 0, one obtains 6(®,, 4)/ô(r, 0) = 0, 


and hence ®, = ®,(A). Faraday's law now implies E = — (4) VA, leading to 
the equation 


V,-5B, = (rsin6) 9A), (6.4) 
where a = w’(A)/. Since y is solely a function of A (6.1) and (6.2) yield 
=-B (6.5) 


(cf. (3.6)), so that the poloidal components of V and B are parallel. The ¢ 
component of the momentum equation integrates to give 


(r sin 0) y, Ee = L(A), (6.6) 


which is the angular-momentum integral. The energy equation in the fixed 
reference frame (2.9) leads to the integral 


B, D(A) 
pyp P pop 74 80. (A), 6.7 
* (v-0e ° pa | Pn 
whereas the entropy equation (2.3) integrates to give 
p = K(A) pr. (6.8) 


Equations (6.3)-(6.8) are the basic integrals governing the flow, where ®,(A), 
V(A), L(A), E(A) and K(A) are arbitrary functions of A. 

By combining the energy integral (6.7) and the angular-momentum integral 
(6.6), one obtains the equation 


E,o(A) = HV} + (Vs—Q,rsin 6)?]+—"> —0,—10:?sin*ó, — (6.9) 


(y- Up 
where 
QA) = (A) (6.10a) 
defines the angular velocity of the plasma and 
E, (4) = E(A)—9,(A) L(A) (6.105) 


is the energy integral in the rotating frame (see e.g. Sakurai 1985). The right- 
hand side of (6.9) represents the sum of the kinetic energy, enthalpy, 
gravitational and centrifugal potential energy of the plasma. Note that in the 
rotating frame V is parallel to B. 

It is now straightforward to solve (6.4) and (6.6) for V, and B,, to obtain 


AM*.(L — Q, ? sin? 0) 


V, 2 O,rsinO-- (5 —i)rsinó 


(6.11) 
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ua (L — Q, r” sin? 0) 
where M*, = ua? /p = V2/V,, defines the poloidal Alfvén Mach number. Note 
that for non-singular behaviour at the Alfvén surface (defined by M, = 1), it is 
necessary that (rsin 6), = 7, = (L/Q,)!, where v, is the Alfvén radius from the 
axis of rotation. 
By taking the component of the momentum equation (2.2) parallel to VA 
(note that B. V4 = 0), we obtain the trans-field momentum equation, or 
generalized Grad-Shafranov equation, in the form 


a [ðf e A) 1O( a ôA 
pr’ sin O|dr\psin@ Or] 1? d6\psin 6 00 


1 0/1 04V 10/1 0A 
up” sala ama ar) T] =Q, (6.13) 


where 


pr B a’ Ma (L-9, r? sin? 6) 


Q- BAE a (1—M3yrsnó 


LL'  (L'-Q,r*sin*0) (L—O,r*sin? 0) 


"Pswo' ——Q-M)sswoó —— c 019 
(see e.g. Sakurai 1985; Heyvaerts & Norman 1989). 
6.2. The potential transonic equation 
For the special class of flows for which 
ELÉK-a-L-0,-0 (6.15) 
the trans-field equation (6.13) may be written in the form 
V x [(1—4/3) B,] = 0. (6.16) 


It is now clear from the development of $4 that Gauss’ law V.B, = 0 may be 
written in the form 
V. (ob) = 0, (6.17) 
where 
= B, = 1—M5, 
bm SB diu 1-4} 


and the trans-field equation (6.16) takes the form 
Vxb=0 (6.19) 


characteristic of irrotational compressible fluid dynamics. 

The Bernoulli integral (6.9), as in §4, may be used to define the sound speed 

a of the corresponding pseudo-fluid. For flows satisfying the conditions (6.15), 
differentiation of the energy integral (6.9) leads to the equation 

B k ldo db dó, dw 

0= dE rot (A) = Worte X Es ) ; 


g 


(6.18) 


(6.20) 
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where w = rsin6, (6.21) 
TE Ta (6.22) 

= - = LET ES EU NAT OF (6.23) 

Mig = yË, Moe m=. (6.24) 


In (6.18)-(6.24), M = b/a is the Mach number of the pseudo-fluid with sound 
speed a, w = rsin ĝ is the radial distance from the axis of rotation (the z axis) 
and Vi, = B,/ (up) defines the ¢ component of the Alfvén velocity. The 
formulae (4.14) also apply in the present case, and the pseudo-fluid is supersonic 
in the sense M > 1 in the regions V, > V, and V, « V, < V,, but subsonic in the 
regions 0 < V, < V, and V, < V, < V, where V,, V, and V, are defined in (3.23) 
and (3.24). 
From (6.20), one finds 


10e a[i, ab, ,, ob), 1 00, 
iaz TA ab, . SEE 
oo ee 5 2o v eee 1 
Qu AM EG e zx) | (6:0) 


for the components of Va/c. Using the results (6.25) and (6.26) for Va/o in the 
continuity equation V.(cb) = 0 for b yields the equation 


ðb ôb 0b, ôb \ xb?-a*0b, b? 
2 h2) 2. þh?) —— ptt a So = 
(a? —b?) as * b?) E? ha (e n v o yz V9, 0. 
(6.27) 
The condition V x b = 0, (6.19), implies that b„ and b, satisfy the equation 
Ob, ôb, 
LE £f. 2 
ôz Ow R MO) 
Alternatively, writing b = V¢ results in the equation 
M cse ao b+- |V -+ 
(a — 5) Fant (0 - 99 F228, Fo Feo 9 $2 vg. va, = o. 
p 


(6.29) 


Equation (6.29) is the usual potential transonic flow equation for compressible 
irrotational non-magnetized fluid mechanics (Courant & Friedrichs 1948, chap. 
1, p. 28), supplemented by the gravitational-force term (oc V®,), and rotational 
and magnetic effects (oc y). For M? > 1 the transonic potential equation (6.29) 
is hyperbolic, but for M? < 1 it is elliptic. 


7. Concluding remarks 

In this paper we have explored the relationship between the generalized 
Grad-Shafranov equation or transverse force balance equation for steady MHD 
flows with one ignorable spatial co-ordinate and the potential transonic flow 
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equation. The characteristics of both equations are related to the group velocity 
surface and Mach cone for linear magnetosonic waves emitted from a given 
point in the flow, in which the edge of the Mach cone corresponds to the 
characteristic directions. Imai’s (1960) derivation of the potential transonic 
flow equation for magnetized flows assumed that the magnetic field was aligned 
with the flow. The derivation of the potential transonic flow equation in the 
present analysis follows from setting the effective source term in the generalized 
Grad-Shafranov equation equal to zero, and applies for flows that are not 
necessarily field-aligned. 

Examples of solutions of the MHD potential transonic flow equation with an 
ignorable co-ordinate z for a rectangular Cartesian co-ordinate geometry were 
obtained by means of the hodograph transformation. Solutions corresponding 
to purely radial flow in the (x,y) plane, helical flows and spiral flows were 
investigated. The solutions obtained differ from a similar class of solutions 
investigated by Seebass (1961) in that the present solutions allow for 
components of the magnetic field and fluid velocity along the z axis not included 
in Seebass' analysis. The energy integral for these flows is conveniently 
represented as a potential-well equation in which the potential U(p) depends 
only on the density. For the class of flows considered by Seebass the potential 
energy U(p) equals the gas enthalpy and there is only one section to the well. 
By including the effect of V, and B,, the well splits into two parts, namely a sub- 
Alfvénie well with M} « 1 and a super-Alfvénie well with AM?» 1. Both 
hyperbolic flows in which the effective Mach number M > 1 and elliptic flows for 
which M? < 1 are obtained. For hyperbolic flows the Mach number M = V,/c, 
where V, is the flow speed in the (x,y) plane and c is the appropriate 
magnetosonie speed on the Mach cone (the interpretation of M as a Mach 
number for elliptic flows is more problematic, since in some flows M? « 0). The 
Mach number M based on the generalized Grad-Shafranov equation charac- 
teristies is also the Mach number of the pseudo-fluid in Imai's analysis. 

There are in general two possible super-Alfvénic flows that join at an inner 
limiting cylindrical radius r = r, which corresponds to a limit line of the 
hodograph transformation. The flows are restricted to the region r > r,. The 
limit line occurs at the point in the flow where the radial flow speed V, matches 
the phase speed c of magnetosonic waves on the Mach cone. One branch starting 
at r = r, consists of a superfast magnetosonic flow (V, > Vj) that accelerates 
with increasing r in which the density p—>0 as r— oo. The other solution 
starting from r = r, consists of a near-Alfvénic decelerating flow solution with 
M , > 1. Similarly, in the sub-Alfvénic flow regime there are two possible flows 
outside the limit circle r = r,,. The two sub-Alfvénic flow solutions consist of a 
hyperbolic flow regime at small radii in which the flow speed V, lies in the range 
V, « V, « V, (where V, is the slow magnetosonie wave phase speed for waves 
propagating along the projected magnetic field B, in the (x, y) plane; if B, = 0, 
V, corresponds to the slow magnetosonic group velocity for waves propagating 
perpendicular to B) followed by an elliptic flow regime at larger r. 

Superimposing flow potentials for purely radial and circular flow results in 
spiral flow solutions. The near-Alfvénic flow solutions with M, z 1 are found to 
have more pronounced spiral structure than both the high-Alfvénic-Mach- 
number branch solutions (M, > 1) and the low-Alfvénic-Mach-number branch 
solutions (M, < 1). Both the projections of the streamlines and the charac- 
teristics in the (x,y) plane and in the (b,,b,) hodograph plane have been 
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investigated. In the hodograph plane the super-Alfvénic characteristics for 
hyperbolic flow are similar to the epicycloid characteristics of pure gasdynamics 
(see e.g. Courant & Friedrichs 1948, chap. 4). However, the characteristics are 
distinctly different than both those of pure gasdynamics and those for MHD 
flows with V, = B, = 0 (Jeffrey & Taniuti 1964), in that a section of the 
characteristics extends inside the limit circle. Other solutions of the potential 
transonic flow equation are also possible, including magnetized versions of 
Ringleb flows and simple waves. 

It is clear that the potential transonic flow equation with an ignorable co- 
ordinate can be obtained for other geometries. In particular, a derivation of the 
potential transonic flow equation and its relation to the generalized Grad—- 
Shafranov equation for the case of flows with an ignorable azimuthal co- 
ordinate ¢ of a spherical polar co-ordinate system (r, 0, 9) has been carried out 
in $6. This geometry is probably the most suitable for models of magnetized 
winds (see e.g. Heinemann & Olbert 1978; Sakurai 1990). 

Solutions of the generalized Grad-Shafranov equation are useful for the 
elucidation of the propulsion and collimation of winds and jets (see e.g. 
Camenzind 1986, 1987; Heyvaerts & Norman 1989; Sakurai 1990). Solutions of 
the potential transonic flow equation are useful for explaining aspects of MHD 
flow past planetary bodies (see e.g. Spreiter et al. 1966, 1970; Spreiter & Rizzi 
1974), sub-Alfvénic flows in the solar atmosphere (see e.g. deVille & Priest 1991) 
and aspects of wind flows and jets. Related work by Sonnerup & Hau (1994) on 
field-aligned flows in CGL plasmas and by Rutkevich & Mond (1992) on 
magnetosonic wave propagation in spiral magnetized flows including shocks are 
other areas of application. 
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Appendix 
The coefficient d in the trans-field equation (3.20) takes the form 


V2 
d= EA) y R- [VS — Và) -A] 


p Ap 


EE 2 a1 m Vx 24,7 V3) 


a Vp— Vis V 
, HPV | ` TELLE | l liea) 
+A Ael VV? Vin | Oot eA +A) Ate 
vive 20 A(A—0%) VÀ 
+O; up V: EE —Vi(A—05)——  —É 42 2 
PAP em (AT Mp) RY 


K'(A)  upVioc A 7 E 
"roO TG Va» ny V:VA.V6,. (A 1) 
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